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Abstract 

We present an exact collapsing solution to 2+1 gravity with a nega- 
tive cosmological constant minimally coupled to a massless scalar field, 
which exhibits physical properties making it a candidate critical solu- 
tion. We discuss its global causal structure and its symmetries in re- 
lation with those of the corresponding continously self-similar solution 
derived in the A — case. Linear perturbations on this background 
lead to approximate black hole solutions. The critical exponent is 
found to be 7 = 2/5. 
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1 Introduction 



Since its discovery, the BTZ black hole solution [1] of 2+1 dimensional AdS 
gravity has attracted much interest because it represents a simplified context 
in which to study the classical and quantum properties of black holes. A 
line of approach which has been opened only recently 0, ||, ||] concerns 
black hole formation through collapse of matter configurations coupled to 
2+1 gravity with a negative cosmological constant. As first discovered in 
four dimensions by Choptuik Q, collapsing configurations which lie at the 
threshold of black hole formation exhibit properties, such as universality, 
power-law scaling of the black hole mass, and continuous or discrete self- 
similarity, which are characteristic of critical phenomena 0. In the case of 
a spherically symmetric massless, minimally coupled scalar field, a class of 
analytical continously self-similar (CSS) solutions was first given by Roberts 
|], ^, [U]]. These include critical solutions, lying at the threshold between 
black holes and naked singularities, and characterized by the presence of null 
central singularities. Linear perturbations of these solutions |ll|, [C| lead to 
approximate black hole solutions with a spacelike central singularity. 

Numerical simulations of circularly symmetric scalar field collapse in 
2+1 dimensional AdS spacetime were recently performed by Pretorius and 
Choptuik [|| and Husain and Olivier ||. Both groups observed critical 
collapse, which was determined in [Q] to be continuously self-similar near 
r = 0. In Q, Garfinkle has found a one-parameter family of exact CSS 
solutions of 2+1 gravity without cosmological constant, and argued that 
one of these solutions should give the behaviour of the full critical solution 
(A ^ 0) near the singularity. 

The purpose of this paper is to present a new CSS solution to the field 
equations with A = which can be extended to a threshold solution of 
the full A ^ equations. The new A = solution is derived in Sect. 3. 
It presents a null central singularity and, besides being CSS, possesses four 
Killing vectors. In Sect. 4 we address the extension of this CSS solution to a 
quasi-CSS solution of the full A < problem, and show that the requirement 
of maximal symmetry selects a unique extension. This inherits the null 
central singularity of the A = solution, and has the correct AdS boundary 
at spatial infinity. Finally, we perform in Sect. 5 the linear perturbation 
analysis in this background, find that it does lead to black hole formation, 
and determine the critical exponent. 
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2 CSS solutions 



The Einstein equations for cosmological gravity coupled to a massless scalar 
field in (2+1) dimensions are 

Gfj, v - J^g^u = kT^ u , (2.1) 
with the stress-energy tensor for the scalar field 

Tfj, v = d^d^ - ^ u d x <pd x <p . (2.2) 

The signature of the metric is (H ), and the cosmological constant A is 

negative for AdS spacetime, A = — l~ 2 . Static solutions of these equations 
include the BTZ black hole solutions fl]] with a vanishing scalar field <f> = 0, 
and singular solutions when a non-trivial scalar field is coupled with the 
positive sign for the gravitational constant k [13]. 

We shall use for radial collapse the convenient parametrisation of the 
rotationally symmetric line element in terms of null coordinates (u,v): 

ds 2 = e 2a dudv - r 2 d6 2 , (2.3) 

with metric functions cr(u, v) and r(u, v). The corresponding Einstein equa- 
tions and scalar field equation are 





(2.4) 


2cr w = ^e 2a - K<f>, u <j) jV , 


(2.5) 


2<7 \u r ,u ~ r ,uu = Kr( P%^ 


(2.6) 


,vT ,v r ,vv = Krcf) 2 ^^ 


(2.7) 


r<t>,uv + r tU (j) )V + r tV (j> )U = . 


(2.8) 



From the Einstein equations, the Ricci scalar is 

R = -6A + 4ne- 2a (j) tU (f) tV . (2.9) 



It follows from ( |2.9| ) and (2.5) that the behavior of the solutions near the 
singularity is governed by the equations (|2.4j )- (|2.8| ) with vanishing cosmo- 
logical constant A = (see also Q). Assuming A = 0, Garfinkle has found 
j|] the following family of exact CSS solutions to these equations 

ds 2 = _ A (W^±VE^t\ dudv-hv + ufde 2 , 

\ -uv J 4 V 
4> = -2cln(^+ v /Z u), (2.10) 
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depending on an arbitrary constant c and a scale A > 0. In ( [2. 10 ), 



u is 

retarded time, and —v is advanced time. These solutions are continuously 
self-similar with homothetic vector (ud u + vd v ) . An equivalent form of these 
CSS solutions, obtained by making the transformation 

- u = (-uf\ v=(v) 2q (l/2q = l- KC 2 ) (2.11) 

to the barred null coordinates (u, v), is 

ds 2 = + (-u)«) 2 ( 2 «- 1 V« dedw - ^(w 35 - (— S) 25 ) 2 ci© 2 , 

<f> = -2c\n{v q + {-u) q ) . (2.12) 

The corresponding Ricci scalar is 

R = ^-{v q + {-u) q f( l - Zq)/q (-u) q - l {v) q - 1 . (2.13) 



Garfinkle suggested that the line element ( [2.10 ) describes critical collapse 



in the sector r = — {u + v)/2 > 0, near the future point singularity r = 
(where the Ricci scalar behaves, for v oc u, as u~ 2 ). The corresponding 
Penrose diagram (Fig. 1) is a triangle bounded by past null infinity u — > — oo, 
the other null side v = 0, and the central regular timelike line r = 0. For 
kc 2 > 1 (q < 0), the Ricci scalar 

i? ~ (y)"-i ~ (^)("- 1 )/ 2 '? (2.14) 

is regular near f = 0, which moreover turns out to be at infinite geodesic 
distance. To show this, we consider the geodesic equation 

( e 2a vJ = -2rr^J 2 = -2/ 2 r" 3 r u (2.15) 

{I constant) near v = 0, u constant, which gives v oc (ls) iq for I ^ 0, or s 2q 
for Z = 0, so that in all cases the affine parameter s — > oo for v — > 0, and the 
spacetime is geodesically complete. For kc 2 < 1 (q > 0), we see from ( p. 13 ) 



that the null line v = is a curvature singularity if kc 2 < 1/2 (q < 1). If 
1/2 < kc 2 < 1 (q > 1), the surface v = is regular. However, as discussed 
by Garfinkle, the metric ( |2.12| ) can be extended through this surface only for 
q = n, where n is a positive integer. For n even, the extended spacetime is 
made of two symmetrical triangles joined along the null side v = 0, and has 
two coordinate singularities r = 0, one timelike (u—v = 0) and one spacelike 
(u + v = 0), but no curvature singularity. For n odd, one of the r = sides 
becomes a future spacelike curvature singularity (e 2cr = 0), similar to that of 
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Brady's supercritical solutions for scalar field collapse in (3+1) dimensions 
Hi, except for the fact that in the present case the singularity is not hidden 
behind a spacelike apparent horizon (Fig. 2). 



Let us point out that, besides the solutions (|2.10[) , the system (|2.4|) -(^8| 



also admits for A = another family of CSS solutions 

1 




<1«- - A\ — - [ du dr ~^(v + u) 2 d9 2 , 

= Ai^-i-uy^-^dudv-^i^-^)^) 2 ^, (2.16) 

with <fi = — 2cln(y / w — \/—u), and we choose A > and consider the sector 
< v < —u. These solutions have a future spacelike central (r = 0) curva- 
ture singularity at (— u) q = v q (where the Ricci scalar ( |2.13| ) diverges) for 
all q < or q > (implying q > 1/2). For q < 0, the Penrose diagram is 
a triangle bounded by past null infinities u — > — oo and v = (which is at 
infinite geodesic distance). For q > 0, geodesies terminate at v = 0, unless 
q = n integer. For n even, the extended spacetime has two central curva- 
ture singularities r = 0, one spacelike and the other timelike. The extended 
spacetime for n odd is more realistic. In this case the extension from v > 
to v < amounts to replacing ( |2.16 ) with A > by the original Garfinkle 



solution ( |2.10| ) with A > 0, the resulting Penrose diagram being that of Fig. 
2. 



3 A new CSS solution for A = 

Among the one-parameter (c or q) family of CSS solutions ( |2.10| ), the special 
solution, corresponding to kc 2 = 1, 

ds i = a(V^+V^1) 4 —- - t(« + uf de 2 , (3.i) 

u v A 

is singled out by the fact that the transformation ( |2.11| ) breaks down for 
this value. The transformation appropriate to this case, 

-u = 2e- u , v = 2e v = 2e u ~ 2T (3.2) 

(with T > U for u + v < 0) transforms the solution (|3.l| ) to 

ds z = e- 2U [-AA(l + e u - T ) 4 dUdV-(l-e 2 ^ u - T ^ 2 de 2 }, 
4> = U- 21n(l + e u ~ T ) (3.3) 
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(we use from now on units such that k = 1, and have dropped an irrelevant 
additive constant from <f>). 

Starting from this special CSS solution of the Garfinkle class, we now 
derive, by a limiting process, a new CSS solution which, as we shall see, 
exhibits a null singularity. We translate T to T — Tq, and take the late-time 
limit To — > — do, leading to the new CSS solution (written for A = —1/2) 

ds 2 = e - 2U {2dUdV - d6 2 ) , 4> = U, (3.4) 

with a very simple form which is reminiscent of the Hayward critical solution 
for scalar field collapse in 3+1 dimensions 12], 

ds 2 = e 2p (2dT 2 - 2dp 2 - dn 2 ) , <f> = r. (3.5) 

The transformation 

u = -e~ 2U , v = V (3.6) 



leads from (3.4) to the even more simple form of this solution 

ds 2 = dudv + udO 2 , (f) = — - ln(— u) , (3-7) 

which is reminiscent of the other form of the Hayward solution 

ds 2 = 2 dudv + uvdn 2 , = -- ln(-u/v) . (3.8) 

The solution (3.4) or (3.7) is continuously self-similar, with homothetic 
vector 

K = d v = -2ud n • (3.9) 
It also has a high degree of symmetry, with 4 Killing vectors 

Lx = du + 2Vdv + 0d e , 
L 2 = 9d v + Ud e , 
L3 = d v , 

U = d , (3.10) 

generating the solvable Lie algebra 

[L 1 ,L 2 ]=L 4 -L 2 , [L 2 ,L 3 ]=0, 
[Li, L3] = — 2L 3 , [L 2 ,L 4 ] = — L 3 , 
[L 1 ,L 4 ] = -L 4 , [L 3 ,L 4 ]=0. (3.11) 
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The Ricci scalar ( |2.9| ) is identically zero for the solution (|3.4|), for which 
the sole nonvanishing Ricci tensor component is Ruu = 1- R follows that 
this metric is devoid of curvature singularity. However there is an obvious 
coordinate singularity at U — > +00, or u = (where r = 0). To determine 
the nature of this singularity, we study geodesic motion in the spacetime 
(|3.7|) , The geodesic equations are integrated by 

^ = vr, u6 = l, TTV + W = e, (3.12) 

where ir and I are the constants of the motion associated with the Killing 
vectors L3 and L4, and the sign of e depends on that of ds 2 along the 
geodesic. The null line u = can be reached only by those geodesies with 
7r / 0. Then, the third equation ( |3.12 ) integrates to 



el el 2 ,, 
v = — ttU 9 + const. = —^u ~ m(— it) + const. . (3.13) 



IT 7T 7T 7T 



It follows that nonradial geodesies (I ^ 0) terminate at u = 0,v — > +00, 
while radial geodesies (Z = 0), which behave as in cylindrical Minkowski 
space, can be continued through the null line u = to u — > +00 . So in 
this sense only the endpoint v — * +00 of the null line u = is singular. 
However formal analytic continuation of the metric ( |3.7| ) from u < to 

u > involves a change of signature from (H ) to (H h), leading to the 

appearance of closed timelike curves. So the null line u = corresponds to a 
singularity in the causal structure of the spacetime, analogous to the central 
singularity in the causal structure of the BTZ black holes . The resulting 
Penrose diagram, reminiscent of that of the Hayward critical solution |l2j , is 
a diamond bound by three lines at null infinity (v = —00, u = —00, v = +00) 
and the null singularity u = (Fig. 3). 



4 Extending the new solution to A 7^ 

In the preceding section we have found an exact solution for scalar field 
collapse with A = 0, which presents a central null singularity. This prop- 
erty makes it a candidate threshold solution, lying at the boundary between 
naked singularities and black holes. However black holes exist only for A < 0, 
so the solution ( |3.7| ) can only represent the behavior of the true threshold 
solution near the central singularity, where the cosmological constant can be 
neglected. This hypothetical A < solution cannot be self-similar, essen- 
tially because the scale is fixed preferentially by the cosmological constant 
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[0. So what we need is to find some other way to extend (3.7) to a solution 
of the full system of Einstein equations with A < 0. 

A first possible approach is to expand this solution in powers of A, with 
the zeroth order given by the CSS solution ( |3.7D . In the parametrisation 
(p. 3D , this zeroth order is (dropping the bars in (13.1 



r = (-u) 1 / 2 , ^ = 0, <h = --hi\u\. (4.1) 

We look for an approximate solution to first order in A of the form 

r = {-u) 1/2 + Ar x , cr = Ao-i, (j> = ~ In |u| + A0i, (4.2) 

with the boundary condition that the fonctions r±, o\ and <j)\ vanish on the 
central singularity u = . Eq (^J|) gives 

r 1 = (- u ) 1 / 2 (±uv + f{u)), (4.3) 
with /(0) = 0. Then, the linearized Eq. fl2.7j ) gives 

2ro /2 (ry 2 0i,J,, = -ri,M,u = \(-u) X/ \ ( 4 - 4 ) 
which is solved by 



{—uv + g{u)). (4.5) 



The linearized Eq. (|2.5|) 



^l,uv = 1 - <Po,u<Pl,v = ( 4 - 6 ) 

then gives 

4 

a 1 = — U v + h(u). (4.7) 
15 



Finally Eq. (2.5) leads to the relation between the arbitrary functions /, g, 
h 

uf"( U ) + f , (u)=g'(u)+ti(u). (4.8) 

Not only does this first order solution break the continuous self-similarity 
generated by ( |3.9|) , as expected, but it also breaks the isometry group gen- 
erated by the Killings (3.10) down to £7(1) (generated by L4 = dg), except 



in the special case / = g = h = 0, where the Killing subalgebra (Li,!,^) 
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remains. This suggests looking for an exact A < extension of the A = 
CSS solution of the form 

ds 2 = e M^dudv + up 2 (x)d9 2 , (j> = --ln\u\+^(x), (4.9) 

with x = uv. This will automatically preserve to all orders the Killing 
subalgebra (Li, L 4 ). Inserting this ansatz into the field equations (|2.4|) - (|2.8[) 
leads to the system 



xp +-p 
2(xa" + cj') +^'{xil>'- 



2a 



x\-p" + 2p'a' - pip' 2 ) + x(-p' + p(o' + tf)) 

-p" + 2pV - pty 



r> 



2x(piP')' + -ftf 



A 

A 

2 L 




v> 



2<T 



(4.10) 

(4.11) 

(4.12) 
(4.13) 

(4.14) 



(' = d/dx). The unique, maximally symmetric extension of the CSS solution 
(3.7) reducing to (|3.7D near u = is the solution of the system (|4.10| )-( 4.14 ) 
with the boundary conditions 



p(0) = 1, (7(0) = o, v(o) = o. 
The comparison of ( |4.12 ) and ( |4.13j ) yields 



p = e CT+ ^. 

The combination ( [4. 1Q| ) + £(4.13) then gives, together with ( 4.16|) , 
x(2a' 2 + 2a'ip' - V' 2 ) + |((/ + i//) = \e 2 ° . 

The third independent equation is for instance ( 4.11| ): 

1, A 



2{xa" + a')+i>\xi>' --) = -e 



2a 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



Using these last two equations with the boundary conditions ( [4.15 ), one can 
in principle write down series expansions for cr(x) and ip(x). Another simple 
relation, deriving from ( |4.13| ) and ( |4.16| ), is 



a" + tf," - cj' 2 + 2iP 



i,l2 



0. 



(4.19) 
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We are interested in the behavior of this extended solution in the sector 
u < 0, v > 0, i.e. x < 0. In this sector, Eqs. (|4.10|) , (f4.14j) and (|4.11[) can 
be integrated to 



-xf 2 p' 



a r° 



(-x) 1/2 pe 2ff dx, 



4 

1 f°,A 



-xcr 



xip'))dx. 



(4.20) 
(4.21) 
(4.22) 



As long as p > 0, Eq. ( |4.20| ) (with x < 0, A < 0) implies // < 0, so that 
p(x) decreases to 1 when x increases to 0. It then follows from ( ^2l| ) that 
i'' < 0. Also, ( PHI) can be integrated by parts to 



xtp' 



16(-x)V4 p 



- 3 / 4 pdx, 



(4.23) 



showing that xip' < 1/4. It then follows from ( 4.22| ) that a' < 0. So, as x 
decreases, the functions p and e 2<J increase and possibly go to infinity for a 
finite value x = x\. If this is the case, the behavior of these functions near 
x\ must be 



p (- + — 

x 4xi 



xln(x) 



48x 2 



+ 



,2<T 



4x 



Ax 2 



x 2 ln(x) 
48xf 



+ 



4xi 



32x 2 



ln(x) + ... 



(4.24) 



(x = X — Xi). 

These expectations are borne out by the actual numerical solution of the 
system 



xp" + -p' 
-p'p + App'a 



p ,2 + p 2 *' 2 



(4.25) 



(this last equation comes from ( 4.13 ) where ip' is given by derivation of 
( |4.16| ) ) where we have set A = —2, with the boundary counditions p(0) = 1, 
p'(0) = -2/3 (see eqs. (fO) and (|J) ), <r(0) = 0. The plots of the functions 
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p(x) , a(x) and ip'(x) are given in Figs. (4,5,6,). The value of x\ is found 
to be approximately —1.94 (i.e. Ax\ = +3.88). 
The coordinate transformation]] 



u 



K- l e~°, v = e V (U = T-R, V = T + R) (4.26) 



leads to x = A 1 e 2R and, on account of ( |4.9[ ) and ( 4.16| ), to the form of the 
metric 

& 2 = - K -^ e ^iR)+m {d tj d y _ e 2m)-v d9 2y {4 2J) 

Near the spacelike boundary R = Ri of the spacetime, the collapsing metric 
and scalar field behave, from ( 4.24| ), as 

ds 2 ~-A -1 (i?i- R)- 2 {df 2 -dR 2 -e Tl - T d9 2 ), </> = 0i + T/2 (4.28) 

(R — R ~ x/2xi). This metric is asymptotically AdS, as may be shown by 
making the further coordinate transformation, 

R-R 1 = -2/XT, T — T\ = 21n(T/2) , (4.29) 

leading to 

ds 2 ~ -A" 1 (^X 2 dT 2 " ^ " X 2 d9 2 ^j , = 0! + ln(r/2) . (4.30) 

The next-to-leading terms in the metric containing logarithms, this asymp- 
totic behavior differs from that of BTZ black holes. 

It follows from this discussion that the Penrose diagram of the A < 
threshold solution in the sector v > 0, u < is a triangle bounded by the 
null line v = 0, the null causal singularity u = 0, and the spacelike AdS 
boundary X — > oo. The null singularity u = remains naked, i.e. is not 
hidden behind a trapping horizon, which would correspond to 

d v r = -(-uf /2 p'{x) = 0, (4.31) 

because p' < (as discussed above) implies that the only solution of this 
equation is u = 0. 

For the sake of completeness, let us also discuss the behavior of the 
solution of the system ( 4.10[ )-( [4.14 ) in the sector x > 0. In this case, one 
can write down integro-differential equations similar to ( |4.20 )-( 4.22 ), from 

x We have taken care that in (4.9) u has the dimension of a length squared while v is 
dimensionless. 
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which one again derives that p' < 0, ip' < and a' < 0. It follows that 
the metric function e 2a decreases as x increases, eventually vanishing for a 
finite value x = xq, corresponding to a spacelike curvature singularity (this 
has been confirmed numerically). The behavior of the solution near this 
singularity is found to be 

7 2 r 
tp ~ 7ln(x - x), a ~ — ln(x - x), p oc (x - x) (7 = v3 — 1), 

(4.32) 

and the coordinate transformation u = e u ,v = e v (x = e 2T ) leads to the 
form of the metric near the singularity 

ds 2 ~ (T - Ty\dT 2 - dR 2 ) + e Ro ~ R (T - T) 2 d6 2 . (4.33) 



5 Perturbations 

To check whether the quasi-CSS solution (^1]) of the full A / problem 
determined in the preceding section is indeed a threshold solution, we now 
study linear perturbations of this solution. Our treatment will follow the 
analysis of perturbations of critical solutions in the case of scalar field col- 
lapse in 3+1 dimensions [11, 12]. 

The relevant time parameter in critical collapse being the retarded time 
U = — (1/2) ln(— u) (the "scaling variable" of |TI[] ), we expand these per- 
turbations in modes proportional to e kU = (—u)~ k ^ 2 , with k a complex 
constant. We recall that only the modes with Re k > grow as U — > +oo 
(u — > —0) and lead to black hole formation, whereas those with Re k < 
decay and are irrelevant. The other relevant variable is the "spatial" coor- 
dinate x = uv , and the perturbations are decomposed as 

r = (- u )V 2 {p(x) + {-u)- k ' 2 r(x)), 

4> = -hn\u\+iP(x) + (-u)- k / 2 4>(x), (5.1) 

a = a(x) + {-i±y k/2 d(x). 

Then, the Einstein equations ( |2.4[ )-(^8|) are linearized in f, <fi, a, using 

5<t>,u = -(-ur k/2 -\x4>i - 5^ v = -(- u )- fc / 2 +V . (5.2) 



The resulting equations are homogeneous in u, which drops out, and the 
linearized system reduces to 

xr" + (-k/2 + 3/2)r = ^e 2<T (f + 2per), (5.3) 
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2xd" + (-jfe + 2)a' = Ae 2a a - (2xifj' - 1/2)0' + (k/2)^'4>, (5.4) 
-(-k + l)xr + ((-k + l)xa' - (k 2 - l)/4)f + pxa' - k(xp + p/2)a = 

-p(x(f>' - k(l/2 - xip')4>) + (1/4 - xfj/)f, (5.5) 
2{p'a' + a'r') - r" = Tp'(2p4>' + ijj'r), (5.6) 
2xp4>" + (2xp' + (-jfe + 5/2)p)ij>f - (k/2)p'4> + (2xip' - l/2)f' 

+ (2xV" + (-fe/ 2 + 5/2)^)f = 0. (5.7) 

What is the number of the independent constants for this system? The 
perturbed Klein-Gordon equation fl5.7p is clearly redundant, while Eqs. (|5.5|) 
and ( p\^ ) are constraints. So, as in the (3+l)-dimensional case [|ll|, ^2]| , 
the order of the system is four, and the general solution depends on four 
integration constants. However, one of these four independent solutions 
corresponds to a gauge mode and is irrelevant. The parametrisation ( |4.S| ) is 
invariant under infinitesimal coordinate transformations v — » v + f(v). For 
f(v) = — av 1+fc / 2 ) these lead to x — » x — a(— ti) _fe / 2 (— x) 1+fc//2 , giving rise to 
the gauge mode 

f fc (x) = a(-x) 1+fc /y(aO , 

4(x) = a(-x) 1+fc /V(x) , (5.8) 
a fc (x) = a[(-x) 1+fc /V'(x)-^(-x) fc /2], 

which solves identically the system (|5.3|)-( |577D . So, up to gauge transforma- 
tions, the general solution of this system depends only on three independent 
constants. 

These will be determined, together with the possible values of k (the 
eigenfrequencies) by enforcing appropriate and reasonable boundary condi- 



tions. We shall use here the "weak boundary conditions" of [12| on the 
boundaries u = and x = x\ (X — > oo) 

lim r _1 ^ 0, lim r + 0, (5.9) 
together with the condition 

f(0) = 0, (5.10) 

which guarantees that the singularity of the perturbed solution starts smoothly 
from that of the unperturbed one. On the third boundary v = 0, we shall 
impose a stronger condition by requiring that the perturbations are analytic 
in v, in order for the perturbed solution to be extendible beyond v = to 
negative values of v at finite u. 
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First, we consider the region x — » where, according to Eqs. (4.1), ( |4,3| ), 
(PDandQ, 

p ~ 1 + -Ax , e 2ff ~ 1 + —Ax , V - —Ax . (5.11) 
3 15 15 

Let us assume a power-law behavior 

f(x) ~ a(-x) p (5.12) 

where p is a constant to be determined. Then Eqs. (|5.3|), fl5,4j ) and (|5.6| ) 
can be approximated near x = as 

xf" + {-k/2 + 3/2)f' ~ Act, (5.13) 
xct" + (-A;/2 + 1)ct'~ -0' (5.14) 
2p'a' ~2pip'4>'. (5.15) 

Eliminating the functions <r and </> between these three equations and using 
Eq. ( 5.11 ), we obtain the fourth-order equation 

4xV" + (-4k + 13)xr " + {k/2 - l)(2Jfe - 5)f" ~ 0, (5.16) 

which implies the power-law behavior ( |5.12| ) with the exponent p constrained 
by 

Pip - l)(p - - 3/4) (p - fc/2 - 1) = . (5.17) 



Obviously the root p = k/2 + 1 corresponds to the gauge mode (|5.q ) and 
must be discarded as irrelevant. As a consequence the general solution near 
x = can be given in terms of three independent constants as 

fix) ~A + B(-x) + AC(-x) 3 / 4+fe / 2 , (5.18) 
a(x) ~ -± + A" 1 - + ^(-x)- 1 / 4 ^ 2 , (5.19) 

0(x) ~ - A" 1 + 5C (fc + 3 )( _ xrl/4+fc/2(5>20) 

Let us note that this solution remains valid in the limit A —* 0, leading to 
the limiting solution f ~ A + B{—x) (with 1? = for k ^ 3), which could 
also be obtained directly by solving the equation f" = which results from 
( |5.6| ) in the limit A — » 0, together with the stronger condition (from Eq. 
(PD) (fc-3)f = 0. 
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Now we enforce the boundary conditions at x = 0. For k > 0, f is 
dominated by its first constant term in ( 5.18| ), so that the condition ( [5 . 1 Q| ) 
can only be satisfied for u — ► if 



A = 0. 



(5.21) 



Then, for k > 1/2, f is dominated by its second term —Bx, leading to a 
perturbation (— u) 1 ^ 2 ~ k ^ 2 f(x) which blows up as u — » and violates ( |5.9| ) 
unless 

A; < 3. (5.22) 

Then we impose the condition of analyticity in v at fixed u. This is satisfied 
if 

fc = 2n-3/2, (5.23) 

where n is a positive integer. Combining eqs. (5.22) and (5.23) we find that 
k has only two positive eigenvalues 



k = l/2, k = 5/2. 



(5.24) 



However, in the above analysis we have disregarded the fact that k = 1/2 
is a double root of the secular equation ( 5.17] ). For k = 1/2 the correct 
behavior of the general solution near x = is 



r(x) ~ A + B(-x) + AC(-x) In \x 
a(x) r 



A A _i5£ 9C 

. . A A ,55 9C 5C, 
(x ) ~ - + A" 1 — + — + — In x , 
4 4 4 4 



(5.25) 
(5.26) 

(5.27) 



which satisfies the condition of analyticity only if C = 0. 

At the AdS boundary (x — * x\) the leading behaviour of the background 
is, from Eqs. fl4.24| ), 



P 



Pi 



e 2 ° ~ 



4rci 



x — x\ \ A J (x — xi) c " 

We again assume a power-law behavior 

or ~ bx q 



(5.S 



(5.29) 



(x = x — xi). Then Eq. (5^), where can be neglected, gives 

</(«-!) = 2, (5.30) 
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i.e. q = — 1 or q = 2. Then, Eq. Q5.3| ) reduces near x = to 

r" - 2x~ 2 ? ~ 4b Pl x q - 3 . (5.31) 

If q = —1, the behavior of the solution is governed by the right-hand side, 
i.e. f oc x~ 2 , which violates the boundary condition (|5.9| ) for x — > X\. So 
the behavior a ~ 6x _1 must be excluded, which fixes another integration 
constant D = (where D is a linear combination of B and C). Then, the 



generic behavior of the solution of Eq. (5.31) with q = 2 is governed by that 
for the homogeneous equation, i.e. 

E 

r . 5.32 

x — X\ 



This is consistent with the boundary condition ( |5.9[ ), and is an admissible 
small perturbation if its amplitude is small enough, E <C p\. 

For k = 1/2, we have seen that two of the three integration constants in 
(|5T25|) -(|5T2l) are fixed (A = C = 0) by condition (|5TC| ) and the analyticity 
condition, while the weak boundary condition at the AdS boundary fixes 
a third constant D = 0. However this is impossible, as the perturbation 
amplitude must remain as a free parameter. So the mode k = 1/2 cannot 
satisfy all our boundary conditions, and we are left with a single eigenmode, 

fc = 5/2, (5.33) 

completely determined up to an arbitrary amplitude by the two conditions 
A = D = 0. 

The corresponding perturbed metric function r behaves near x = as 

r ~ (— u) 1 l 2 [\ + \tix-(-u)-*l*Bx\. (5.34) 
3 

For B < 0, the central singularity r = is approximately given by 

{-u) l/A ~ -Bv. (5.35) 

Our boundary conditions guarantee that it starts at u = v = (as for 
the unperturbed solution) and then becomes spacelike in the region v > 0. 



This singularity is hidden behind a trapping horizon (defined by Eq. ( 4.31 )) 
which, near x = 0, is null, 

(-n) 5/4 = f (5.36) 
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(a null trapping horizon was also found in |12[| ). Let us point out the crucial 
role played by the cosmological constant A in the formation of this trapping 
horizon. For A = 0, p(x) = 1, while, as discussed after Eq. ( [5.20 ), the 



perturbation f with the boundary condition ( |5.10| ) vanishes for A = 0, so 
that the perturbed radial function r is (as in [Q]) identical to the CSS one, 
and the trapping horizon does not exist. Near the AdS boundary x — > x%, 
it follows from ( 5.28 ) and ( 5.32; ) that both the central singularity and the 



trapping horizon are tangent to the null line 

(-) 5/4 = "£(^r 1/2 • (5-37) 

Thus, perturbations of the quasi-CSS solution lead to black hole for- 
mation, showing that this solution is indeed a threshold solution, and is a 
candidate to describe critical collapse. Near-critical collapse is characterized 
by a critical exponent 7, defined by the scaling relation Q oc \p — p*\ sl , for a 
quantity Q with dimension s depending on a parameter p (with p = p* for 
the critical solution). Choosing for Q the radius tah of the apparent hori- 
zon, and identifying p — p* with the perturbation amplitude B, we obtain 
from flOp 

tab - (^) 2/5 , (5.38) 
leading to the value of the critical exponent 7 = 2/5, in agreement with the 



renormalization group argument [14| leading to 7 = 1/k. 



6 Conclusion 

We have discussed in detail the causal structure of the Garfinkle CSS so- 



lutions ( 2.10 ) to the A = Einstein-scalar field equations. From a special 
solution of this class, we have derived by a limiting process a new CSS solu- 
tion, which we have extended to a unique solution of the full A < equations, 
describing collapse of the scalar field onto a null central singularity. This is 
not a curvature singularity (all the curvature invariants remain finite), but 
a singularity in the causal structure similar to that of the BTZ black hole. 
Finally, we have analyzed linear perturbations of the A < solution, found 
a single eigenmode k = 5/2, checked that this mode does indeed give rise to 
black holes, and determined the critical exponent 7 = 2/5. 

For comparison, Choptuik and Pretorius ||J derived, by analysing the 
observed scaling behavior of the maximum scalar curvature, the value 1.15 < 
7 < 1.25 for the critical exponent. This value is different from the value 
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7 ~ 0.81 obtained in the numerical analysis of Husain and Olivier ||] from 
the scaling behavior of the apparent horizon radius. Our value 7 = 0.4, while 
significantly smaller than these two conflicting estimates, is of the order of 
the theoretical value 7 = 1/2 derived either from the analysis of dust-ring 
collapse [15], of black hole formation from point particle collisions 16], or 
of the J = to J 7^ transition of the BTZ black hole flTj . 

It is worth mentioning here that, even though they were obtained for 
a vanishing cosmological constant and thus solve the A / equations only 
near the singularity, the Garfinkle CSS solutions are, for the particular value 
(chosen in order to better fit the numerical curves) c = (7/8) 1 / 2 ~ 0.935, m 
good agreement || with the numerical results of at an intermediate time. 
The fact that this value is close to 1 suggests that the c = 1 CSS solution 
( |3.3[) approximately describes near-critical collapse at intermediate times. If 
this the case, then it would not be surprising if its late-time limit, our new 
CSS solution Eq. ( |3.4D , gives a good description of exactly critical collapse 
near the singularity. A fuller understanding of the relationship between the 
numerically observed near-critical collapse and these various A = CSS 
solutions could be achieved by extending them to A < 0, as done in the 
present work for the special solution (f 
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Figure 3: Penrose diagram of our new CSS solution 
u = is a singularity in the causal structure. 
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Figure 4: Numerical plot of the function p(x) as derived from the system 
(14.250 with p(0) = and p'(0) = -2/3, showing the diver gence of p for 
the AdS boundary is approached (the behaviour is given in the 



first of Eqs. fl424|) ). 
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Figure 5: Numerical graph of a(x) starting from <r(0) = 0. In the limit 
x — > x\ this is well represented in the second of Eqs. ( |4,24| ). 




Figure 6: Plot of ip'(x). In particular it is clear that ip"{x) — > oo as x — > x%. 
This feature is reproduced in the third of Eqs. ( 4.24[ ) (giving ip" ~ ln(x — 

X!)). 
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